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Abstract 

We present a new systematic method of constructing rational mappings as 
ergordic transformations with nonuniform invariant measures on the unit in- 
terval / = [0,1]. As a result, we obtain a two-parameter family of rational 
mappings that have a special property in that their invariant measures can 
be explicitly written in terms of algebraic functions of parameters and a dy- 
namical variable. Furthermore, it is shown here that this family is the most 
generalized class of rational mappings possessing the property of exactly solv- 
able chaos on i", including the Ulam=Neumann map y = Ax(l — x). Based on 
the present method, we can produce a series of rational mappings resembling 
the asymmetric shape of the experimentally obtained first return maps of 
the Beloussof-Zhabotinski chemical reaction, and we can match some rational 
functions with other experimentally obtained first return maps in a system- 
atic manner. 
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Characterizing invariant measures for explicit nonlinear dynamical systems is a funda- 
mental problem which connects dynamical theory with statistics and statistical mechanics. 
In some cases, it would be desirable to to characterize ergordic invariant measures for sim- 
ple chaotic dynamical systems. However, in the cases of chaotic dynamical systems, such 
attempts to obtain explicit invariant measures have rarely been made. One well-known ex- 
ception is the logistic map Y = 4X(1 — X) = fo(X) on I = [0, 1] given by Ulam and von 
Neumann in the late 1940's Jl|. 

The Ulam=Neumann dynamical system Xi = fo(xi-i) has an ergordic measure fi(dx) = 



dx 



■K-^ x(l — X) 



such that the time averages of a function Q(x) can be explicitly computed by the 



formula liniAr^oo J^fLn 1 jiQ( x i) = In Q ^ dx for almost all initial conditions Xn G i\ The 

JV nWx(l— x) 

first attempt to generalize the Ulam=Neumann map within a set of rational functions was 
made by Katsura and Fukuda in 1985 0]. The Katsura and Fukuda model is written as 

for < / < 1. Clearly, the Ulam=Neumann map can be regarded as a special case of 
Katsura=Fukuda systems where the parameter I is set to 0. The author showed || that the 
Katsura=Fukuda mappings ([[]) also have ergordic measures which can be written explicitly 

as 

ii(dx) = p(x)dx = . , (2) 

2K{l)y/x{l-x){l-lx) 

where K(l) is the elliptic integral of the first kind (g = 1) given by K(l) = L , du 

•y/(l— u 2 )(l— lu 2 ) 

It is known f2| that the Katsura=Fukuda systems and the Ulam=Neumann system also have 
explicit solutions in terms of the Jacobi sn elliptic function, as 

x n = sn 2 (K(l)2 n e ), B el, (3) 

where Vl corresponds to the modulus of Jacobi elliptic functions. The validity of the for- 
mulae of the general solutions (|3|) is easily checked using the duplication formula || of the 
Jacobi sn elliptic function 



2sn(u)\/(l — sn 2 (u))(l — lsn 2 (u)) 

^ - V \i-Jl»)) -■ (4) 

Because the Ulam= Neumann system and Katsura=Fukuda systems have not only exact 
solutions (|3]) but also explicitly written ergordic invariant measures (0), we call a dynamical 
system exactly solvable chaos if it has both of these properties. 

Thus, it is of great interest to investigate whether we can generalize the Ulam=Neumann 
system and the Katsura=Fukuda systems further within a set of rational functions, main- 
taining the property of exactly solvable chaos. The main purpose of the present letter is 
to show that by using the addition formulas of elliptic functions we can construct a two- 
parameter family of rational mappings of exactly solvable chaos, and at the same time, that 
there is a certain limitation to generalizing this family within a set of rational functions. 

Our results reported here || concern the following rational transformations 

where A = -2(l+m+lm),B = Sim, C = I 2 + m 2 -2lm-2l 2 m-2lm 2 + l 2 m 2 , and X e I. The 
parameters I and m are arbitrary real numbers satisfying the condition — oo < m < I < 1. 
Figure 1(a) shows various shapes of the proposed mappings @. Surprisingly, some rational 
maps in ([|) strongly resemble the asymmetric shape of the experimentally obtained first 
return maps of the Beloussof-Zhabotinski chemical reaction. Here, we will prove the following 
statement. The two-parameter family of rational mappings (||) is also exactly solvable 
chaos such that the dynamical systems Xi+\ = fi, m {xi) (|5|) have ergordic invariant measures 
explicitly given by 

fi(dx) = p(x)dx = , (6) 

2K(l, m) Jx(l — x)(l — lx)(l — mx) 

where K is given by the integrals 

K(l, m) = f 1 = (7) 

^(l-u 2 )(l-/u 2 )(l-mu 2 ) 

and it has general solutions explicitly given by 



x n = s 2 (K(l,m)2 n 9 ) 6 El. (8) 

We prove this by explicitly computing the duplication formula of the following degenerated 
hyperelliptic function s(x) defined by 

S -U x \ = f X ^ (9) 

^0 ^(1 - y?){l - lv?){l - mv?) 

Here, degenerated means that this hyperelliptic integral of the R.H.S. of Eq. (||]) can be 
reduced to a certain elliptic integral by a rational change of variables. Although the reduction 
of Abelian integrals of genus g > 2 to elliptic functions was intensively studied in the 19th 
century by Jacobi, Weierstrass, Konigsberger, Kovalevskaya and others, it was only in the 
1980's that the theory of the reduction was successfully applied to physics for obtaining, 
for example, explicit periodic solutions in terms of elliptic functions for the Korteweg-de 
Vries (KdV) equation and the sin-Gordon equation ||. Let us consider reduction of the 
hyperelliptic integrals @ as 

s-\ X ) = r du = f dv (io) 

Jo -« 2 )(l-lM 2 )(l-m« 2 ) Jo 2y/v(l - v)(l - lv){\ - m«)' 

where u 2 = v. Thus, we can write the degenerated hyperelliptic function s(x) in terms of 
the Weierstrass elliptic functions. The Weierstrass elliptic function p(u) of u G C is defined 
by 

1 ' 1 1 

= ^ + 5 { (u-2j^-2A^ 2 ) 2 " (2^ 1 + 2^ 2 ) 2} ' (U) 

where the symbol J2' means that the summation is made over all combinations of integers j 
and k, except for the combination j = k = 0, and 2ui and 2u2 are periods of the function 
p{u) 0]. The Weierstrass elliptic function p(u) satisfies the differential equation 

(^) 2 = - g 2 p(x) ~ 9s, (12) 

with the invariants 34(^1,^2) = 60Ej-, fc (M+ 1 fcu;2) 4 and gz{ux,u 2 ) = MOE^ ^j^ 1- Let 
e±, e 2 and be the roots of the equation 4z 3 — g 2 z — g 3 = 0; that is, 



ei + e 2 + e 3 = 0, e x e 2 + e 2 e 3 + e 3 ei = eie 2 e 3 = ^. (13) 

The discriminant A of the function jp(w) is given by A = g\ — 21 g\. If A > 0, all roots e±, e 2 
and e 3 of the equation 4z 3 — g 2 z — g% = are real. Thus, we can assume that e\ > e 2 > e 3 . 
In the case that A > 0, it is known that the periods u\ and uj 2 are written simply as 



/ dz ■ r 3 dz riA\ 

iOi = = , uj 2 = i I . (14) 

Jei V 4 « - #2^ - #3 •'-oo V#3 + 92Z ~ 4z J 

Using the transformation of the variable as v = ^ \7+2m-swl — + 1, we can rewrite s~ 1 (x) 

y 3 

in Eq. (|I~0|) as 

dv f 3 + i- x i — dy s 

(15) 



o 2^/v(l - v)(l - lv)(l - mv) J^P 1 V 4 ?/ 3 - 92V ~ 93 ' 

i A(l-l+l 2 -m+m 2 -lm) j 4(2-Z-m)(2Z-m-l)(2m-Z-l) 1Ir 4-1. j-l. 4- /( 3 

where g 2 = — 3 i and g% = — - — ^ — - -. We note here that Ay° — 

92V — 93 can be factored as 

. ^ 2 — l — m.. 21 — m — l w 2m — Z — 1. 

4?/ 3 - g 2 y ~ 93 = 4(2/ )(y )(y ). (16) 

We set ei = 2 ~'~ m > e 2 = 2 '~^ , ~ 1 > e 3 = gg =^= i . Thus, using the integral representation 
of the period cui (|14D and the differential equation (|12"D for the Weierstrass elliptic function, 
s(x) can be written explicitly in terms of the Weierstrass elliptic function as 



^) = 1 - , (1 ~° (1 2 ,: 2 m) 3 , r - (17) 

p{u x -x) - ^-■• h 



The function s 2 (x) also has the same periods u>i and u 2 computed using formula ([14]). It is 
noted here that because 

A = gl - 21g\ = 16(1 - /) 2 (1 - m) 2 (/ - m) 2 > (18) 

for — 00 < m < I < 1, the period 2cui is always reczi, while the period 2a; 2 is always ptzre 
imaginary. Using the addition formula, 

and the duplication formula, 



(20) 



for the Weierstrass elliptic function p(u) Q, we finally obtain the explicit duplication formula 

of s(x) as 



s 2 (2x) 



As 2 (x)(l - s 2 (x))(l - ls 2 (x))(l - ms 2 (x)) 



(21) 



1 + As*(x) + Bs e (x) + Cs 8 (x) 

where A = —2(1 + m + Im), B = 8lm and C = I 2 + m 2 — 2lm — 2l 2 m — 2lm 2 + l 2 m 2 . If we 
set X = s 2 (x) and Y = s 2 (2x), we obtain system @ as Y — fi >m (X). 
Using the relations 



• 20) = f hm {s 2 {u x 9)), s 2 ^ ■ (2 - 29)) = /,, m (s 2 M)), 



(22) 



for G [0, 1] and by defining the homeomorphism of [0, 1] into itself given by <fii, m (x) 
^s~ l (-s/x), we derive the tent map f(x) = (f) l>m o / Zj7n o as 



/(a:) = 2a; for x G 



1 

°'2 



, f(x) = 2 — 2x for x G 



5.' 



(23) 



Because this tent map fl23|) is clearly ergordic and preserves the Lebesgue measure, the map 
fi tTn preserves the measures 

dx 



fi(dx) 



■'Lm 



dx 



(24) 



dx 2K(l,m)y/x(l-x)(l-lx)(l-mx) 

This measure (|24D is absolutely continuous with respect to the Lebsegue measure, which 
implies that the Kolmogorov- Sinai entropy h(n) is equivalent to the Lyapunov exponent of 
log 2 from the Pesin identity 0, and that the measure ( p4T) is a physical one in the sense 
that it is the Sinai-Ruelle-Bowen (SRB) measure such that for almost all initial conditions 
xo, the time averages liniAr^oo — X^o 1 ~ x «) reproduce the invariant measure fi(dx) [0]. 

In the same way, we can construct generalized cubic maps f lr ^ from the triplication 
formula s 2 (3x) = //m(s 2 (x)) as 

X(-3 + 4X + £ii AiX^ 2 



Y = ffl(X) 



(25) 



where A±, ■ ■ ■ , A4 and B2, ■ ■ ■ , Bg are polynomial functions in the parameters / and m which 
vanish for I — m — ||. The generalized cubic map //^ has the same invariant measures 
(p^) because the relation f^(x) = (j)^ m o o 0^ holds for the piecewise-linear map 
fW(x) = 3x for < x < i -3z + 2 for § < x < § and 3x - 2 for § < x < 1. If we set 
I — m — 0, this rational mapping is reduced to the cubic map Y = X(3 — AX) 2 as a special 
case of Chebyshev maps obtained by Adler and Rivlin . Thus, we can obtain generalized 
Chebyshev maps as rational functions fj^ from the addition formulas s 2 (px) = fi P J l {s 2 (x)), 
which have the same invariant measures (^4|) . The shapes of generalized cubic maps are 
depicted in Fig. 1(b). Based on the case of Beloussof-Zhabotinski map, we predict here that 
some rational mappings (^5l), such as that shown in Fig. 1(b), can resemble the first return 
maps experimentally constructed from some unknown chaotic phenomena. 

Are there more generalized rational mappings that possess the properties of exactly 
solvable chaos such as (|5|)? The exact solvability of the present rational mappings ([5]) and 
(p5|) is due to the fact that s 2 (x) in Eq. ( |17D is a rational function of the Weierstrass 
elliptic function p(u) having the addition formula and the real period ujx. For an arbitrary 
set of parameters e\, e<i{< ei), and e^{= —e.\ — e^) which determine the Weierstrass elliptic 
function p(u) with the real period, there exists a set of parameters / and m of (|5[) given by 
I — 1 — (ei — t2) and m — 1 — (2ei — 2e2) < i.e., the mapping h : (ei,e2) 1— > (l,m) is 
a bijection. In other words, the family {f^ m } has one-to-one correspondence to the set of 
Weierstrass elliptic functions {p(u)} with real periods. Furthermore, any elliptic function 
w(u) can be expressed in terms of Weierstrassian elliptic functions p(u) and p(u)' with the 
same periods, the expression being rational in p{u) and linear in p'{u) pyiQf. Since it is 
known [[11] from the Weierstrassian theorem that the class of functions h(u) having the 



algebraic addition formulae [|T2| including the duplication formulae is exactly the class of 
algebraic functions of elliptic functions p(u), which of course includes the class of algebraic 
functions of sin functions, we can say that this two-parameter family of dynamical systems 
essentially forms a maximal class in representing exactly solvable chaos induced by one- 
dimensional rational mappings. 



In conclusion, we present a new method of constructing ergordic transformations related 
to rational functions with explicit nonuniform invariant measures using the addition for- 
mulas of elliptic functions. As a result, we systematically generalize the Ulam=Neumann 
(logistic) map and Chebyshev maps into two-parameter families of rational mappings. We 
also showed that constructing the more generalized family of rational mappings possess- 
ing explicit ergordic invariant measures has a certain limitation due to the Weierstrassian 
argument concerning the addition formulas for general algebraic functions. As for the appli- 
cations, all of the constructive results given here for ergordic invariant measures by rational 
mappings can be directly used as nonlinear random number generators for the Monte Carlo 
methods. 

The author would like to acknowledge Dr. K. Iguchi and Dr. A. Bobenko for useful 
discussions. Also, Professor T. Kohda and Professor S. Amari provided helpful suggestions. 
This work was supported by the RIKEN Special Researcher's Program to promote basic 
sciences. 



8 



FIGURES 



FIG. 1(a): (mono). The two-parameter family of exactly solvable chaos mappings (|5|) 
for (7, to) = (0.1, —10), (0.1, —20), ••■ (0.1, —100). The asymmetric shapes of this class of 
mappings are very similar to the first return maps of the Beloussof-Zhabotinski chemical 
reaction. 
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FIG. 1(b): (mono). The two-parameter family of exactly solvable chaos mappings 
for (l,m) = (0.1,-10), (0.1,-20),- • -(0.1,-100). 



10 



REFERENCES 

[1] S.M. Ulam and J. von Neumann, Bull. Am. Math. Soc. 53, 1120 (1947). 

[2] S. Katsura and W. Fukuda, Physica (Amsterdam) 130 A, 597 (1985). 

[3] K. Umeno, "Inferring invariant measures of dynamical systems", IEICE Technical Re- 
port NC96-13 (1996). 

[4] E.T. Whittaker and G.N. Watson, A Course of Modern Analysis (Cambridge University 
Press, 1935). 

[5] This work was first reported as an oral presentation; K. Umeno, "Generalized Ulam- 
Neumann difference equations and statistics" , at the International Workshop of Sym- 
metries and Integrability of Difference Equations //held in Canterbury, England (July, 
1996). 

[6] E.D. Belkolos, A.I. Bobenko, V.B. Matveev and V.Z. Enol'skii, Russ. Math. Surv. 
41(1986)1 and references cited therein. 

[7] J. P. Eckmann and D. Ruelle, Rev. Mod. Phys. 57, 617 (1985). 

[8] Precise forms of A\, ■ ■ ■ , A4 and Bi, ■ • • , B g are given as 



11 



A x = 4(1 + m), A 2 = -6(1 + m + lm),A 3 = 12/m, 
A 4 = I 2 + m 2 - 2lm - 2l 2 m - 21m 2 + / 2 m 2 , 
B 2 = -12(l + m + lm), 

B 3 = 8(1 + m + I 2 + m 2 + l 2 m + lm 2 + 15/m), 

B 4 = 6(5/ 2 + 5m 2 - 26/m - 26/ 2 m - 26/m 2 + 5/ 2 m 2 ), 

B 5 = 2A(-2l 2 - 2m 2 - 21 s - 2m 3 + Aim + 7l 2 m + 71m 2 ) 

+24(4/ 3 m + 4/m 3 + 7/ 2 m 2 - 2/ 3 m 2 - 2/ 2 m 3 ), 

B 6 = 4(4/ 2 + 4m 2 + 4/ 4 + 4m 4 + 171 s + 17m 3 - Sim) 

+4(-17/ 2 m - 17/m 2 - 17/ 3 m - 17/m 3 - 8/ 4 m - 8/m 4 ) 

+4(4/ 2 m 4 + 4/ 4 m 2 - 17/ 3 m 2 - 17/ 2 m 3 + 17/ 3 m 3 - 54/ 2 m 2 ), 

B 7 = 2A(—l 3 — m 3 — I 4 — m 4 + l 2 m + lm 2 — l 3 m — lm 3 ) 

+2A(l A m + lm 4 + 4l 2 m 2 + 4l 3 m 2 + 4l 2 m 3 ) 

+24(l 4 m 2 + l 2 m 4 - l 3 m 3 - l 4 m 3 - l 3 m 4 ), 

B 8 = 3(3/ 4 + 3m 4 + 4/ 3 m + 4/m 3 + 4/ 4 m + 4/m 4 - 14/ 2 m 2 ) 

+3(-4/ 3 m 2 - 4/ 2 m 3 - 4/ 3 m 3 - 14/ 4 m 2 - 14/ 2 m 4 + 4/ 4 m 3 + 4/ 3 m 4 + 3/ 4 m 4 ), 

B 9 = 8(-l 4 m - lm 4 + l 3 m 2 + l 2 m 3 + l 4 m 2 + l 2 m 4 - 2l 3 m 3 + l 4 m 3 + l 3 m 4 - l 4 m 4 ). 

[9] R.L. Adler and T.J. Rivlin, Proc. Am. Math. Soc. 15, 794 (1964). 

[10] Thus, any elliptic function w(u) has an algebraic relation with the Weierstrass elliptic 
function p(u) with the same periods as G(w, p) — 0, G being a rational function of two 
variables. 

[11] A. Hurwitz and R. Courant, Vorlesungen iiber allgemeine Funkitionentheorie und ellip- 
tische Funktionen (Springer- Verlag, Berlin, 1964). 

[12] We say that a function w(u) has an algebraic addition formula if there exists a polyno- 
mial function G of three variables such that the relation G(w(ui + u 2 ), w(ui), w(u 2 )) = 
holds for an arbitrary set of the variables U\ and u 2 . 



12 



1=0.1 




1=0.1 




